Quantum phase transition and sliding Luttinger liquid in coupled t — J chains 
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Using a recently proposed perturbative numerical renormalization-group algorithm, we explore 
the connection between quantum criticality and the emergence of Luttinger liquid physics in t — J 
chains coupled by frustrated interactions. This study is built on an earlier finding that at the 
maximally frustrated point, the ground state of weakly-coupled Heisenberg chains is disordered, 
the transverse exchanges being irrelevant. This result is extended here to transverse couplings up 
to J± — 0.6, and we argue that it may also be valid at the isotropic point. A finite size analysis 
of coupled Heisenberg chains in the vicinity of the maximally frustrated point confirms that the 
transverse spin-spin correlations decay exponentially while the longitudinal ones revert to those of 
decoupled chains. We find that this behavior persists upon moderate hole doping x < 0.75. For 
larger doping, the frustration becomes inactive and the quantum critical point is suppressed. 
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I. INTRODUCTION 

In a number of materials including cuprate high tem- 
perature superconductors 1 , organic 2 and inorganic* quasi 
one-dimensional (ID) conductors, the physics is dom- 
inated by the combined effect of low dimensionality, 
strong electron correlations and competing orders. These 
materials display a phenomenology, in the metallic state 
above the ordered states, which departs from the Fermi 
liquid (FL) picture. This has motivated a search for 
a new paradigm which goes beyond FL theory. Two 
new ideas, electron fractionalization 4 *^ and quantum 
criticality 6 (QC), have been proposed as the possible 
physical effects which may lead to the observed non-FL 
behavior. 

The breakdown of the FL theory is known to occur 
in ID where the low energy physics is described by Lut- 
tinger liquid (LL) theory. The LL is fundamentally differ- 
ent from a FL; for instance, instead of the quasiparticle 
peak displayed by the FL at zero energy in the spectral 
weight function, a pseudogap is observed in the LL. This 
is because in the LL the low energy excitations are col- 
lective spin density and charge density, in contrast to the 
FL where they are still electron-like. An interpretation 
of this phenomenon is that the electron fragments into 
density waves known as spinons and holons which prop- 
agate with different velocities. The LL would thus be 
the natural starting point to search for non-FL behav- 
ior in higher dimension. However, it has been impossi- 
ble to find a LL in dimensions greater than one using 
perturbative methods. This has been tried by starting 
from ID models coupled by a small transverse hopping 7 
or directly from the 2D non-interacting electron gas in 
which the interaction is introduced perturbativel}^. It 
was found in both cases that the FL is the stable fixed 
point. An alternative strategy applied recently is to build 
a LL fixed point in anisotropic 2D systems by coupling 
ID LLs by marginal transverse strong forward density- 
density or current-current couplings 9, 10,1 1 . The new 2D 
fixed point, called sliding LL (SLL), retains the physics 
of the ID LL. It was shown that there is a regime of 



parameters in which Josephson, charge and spin density 
waves and single particle hopping are irrelevant, i.e., a 
domain of stability of the 2D LL. There has not been, to 
the author's knowledge, any detailed analysis as to why 
these two approaches lead to conflicting conclusions. It 
would appear that using simple perturbation theory, one 
cannot go smoothly from a single chain to a SLL. 

A different approach has been to assign the departure 
from a FL behavior to the proximity of a quantum critical 
point (QCP) 6 . Because of the proximity of the QCP the 
adiabatic continuity concept, which is the basis of the FL 
theory, is no longer valid. In the normal state, the eigen- 
states of the system are a mixture of states which evolve 
separately into the ground states existing on the two sides 
of the QCP. Because of this mixture, the single-particles 
excitations of the system are not simple electron-like ex- 
citations as in the FL. 

Considering the pure spin one-half limit of the prob- 
lem, various studies have found that any small exchange 
leads to the onset of long-range order. Laughlini^ has ar- 
gued that spin-1/2 systems have the propensity to order 
and it is only at a QCP that they do not. Yet considering 
that experiments clearly show a departure from the FL, 
one can still retain the electron fractionalization hypothe- 
sis as the driving mechanism of the non-FL but attribute 
it not to LL physics but to the presence of a QCP. This 
suggests that in constructing a model Hamiltonian for 
coupled LLs, one must include competing interactions 
that preclude any order at the QCP. 

The difference between the SLL and QCP physics is 
ultimately two different views of the same reality. If 
they are both correct, they should emerge from a well 
controlled study of an appropriate microscopic Hamilto- 
nian. The two-step density-matrix renormalization group 
(TSDMRG)i2*i^ offers a new approach to study these is- 
sues. Although it is a perturbative approach, it has the 
important non-perturbative property that it can lead to 
an ordered state starting from a disordered stat o 14 ' 15 . 
Thus it has the ability to treat effectively the different 
ground states that may arise in weakly coupled chains. 

In this paper, we apply the TSDMRG to study the 
QPT in the spatially anisotropic t — J — t± — J±_ — Jd 
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model. We show that in this model the concept of LL 
can be unified with that of quantum criticality, namely 
that the model undergoes a QPT at Jd = J±/2 (These 
couplings are defined in the Hamiltonian ( ^) given be- 
low). At the QCP the system is in a state made of nearly 
disconnected chains that can be interpreted as a SLL. We 
first analyze the half-filling case which corresponds to the 
anisotropic J\ — J% model. We show that the dimerization 
found in previous approximate field theory analysesi&ii 
does not occur. We then study the doped case at hole 
densities up to x = 0.4. The picture observed in the 
half- filled case remains valid at moderate dopings. At 
the QCP, the transverse spin-spin correlation function 
decays exponentially as in the half-filled case. But at 
higher dopings, the frustration becomes ineffective and 
the transition is suppressed; the ground state is a spin 
density wave at all couplings. This suggests the existence 
of a SLL to FL crossover at finite temperature. 



II. MODEL AND METHOD 

The anisotropic t — J — t± — J± — Jd model is: 
H = -t ^(ct ;C i+u + h.c.) + J ^2 - 2 n hl n i+i 

i,l Z,Z 

-*_L ^2i4d c hi+i + h - c ) + J±- X] s M s M+i _ -^n^im t i + 

i,l 

+ Jd X Sj.fSi+U+1 + Sj+i ) iSj ) i4- 



a 

Jd X SijS i+1;i+1 + S^+i^j £5,^+1,(2) 



where E\\\ n } is the sum of eigenvalues of the different 
chains, E\u n ] = X)z £ n<> l^llfnl) are the corresponding 
eigenstates, |$||[„]) = |0 ni )|0„ 2 )...|0 nz ,); cj;, c M , and 
are the renormalized matrix elements in the single chain 
basis. They are given by 



(5i,i) ni > m > = {-l) ni {<i> ni \c i>l \4> mi ), 
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(3) 
(4) 
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where rii represents the total number of fermions from 
sites 1 to i — 1. For each chain, operators for all the sites 
are stored in a single matrix 
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where t and J > are the in-chain hopping and ex- 
change parameters; tj_ and Jj_ > are the transverse 
hopping and exchange parameters; Jd > is the diag- 
onal exchange parameter. We start by describing the 
extension to fermion systems of the numerical renormal- 
ization method, which has been applied so far only to 
spin system a 13 ! 14 . There is no significant difficulty in 
going from spins to fermions. One should simply be 
careful to consistently apply fermion anticommutation 
rules. It is crucial to choose an order for the sites in 
the 2D lattice and keep it throughout the derivation of 
the algorithm. The method is a special case of a more 
general matrix perturbation method based on the Kato- 
Bloch expansionSJiSi which was recently introduced by 
the author—. The method has two main steps. In the 
first step, the usual ID DMRG method 22 is applied to 
find a set of low lying eigenvalues e n and eigenfunctions 
\4> n ) °f a single chain. In the second step, the 2D Hamil- 
tonian is then projected onto the basis constructed from 
the tensor product of the |</> n )'s. This projection yields an 
effective one-dimensional Hamiltonian for the 2D lattice, 



H 



1,1 



Since the in-chain degrees of freedom have been inte- 
^ grated out, the interchain couplings are between the 
block matrix operators in Eq.( [7| |HJl which depend 
only on the chain index I. In this matrix notation, the 
i(l effective Hamiltonian is one-dimensional and it is also 
studied by the DMRG method. The only difference with 
a normal ID situation is that the local operators are now 
m x m matrices, where m is the number of states kept 
to describe the single chain. In this study, mostly exact 
diagonalization (ED) instead of DMRG is used in the ID 
part of the algorithm. This has the obvious advantage 
that all the eigenvectors and eigenvalues are known up 
to a well defined accuracy which is set to ICP 6 in the 
Davidson algorithm used to obtain them. The two-step 
DMRG method is variationali^i^, so the results can sys- 
tematically be improved by increasing m. In this study, 
the calculations were done on a Dell 670 dual Xeon (3.4 
GHz) workstation with 4GB of RAM. The maximum m 
we can reach for the calculations to run within a rea- 
sonable time (about three days) is about m = 96 for a 
L X (L + 1) lattice with L = 16. 

In this study, we are mostly concerned with the pa- 
rameter regime where Jd ~ Jl/2. It is in this regime 
that the sign problem in the QMC is the most severe. In 
Ref— , we have shown that the two-step DMRG is at its 
best in the highly frustrated regime where QMC simula- 
tions cannot be done. To illustrate this fact, we show in 
Table([IJ the ground state energy (Eq) as a function of 
m for the unfrustrated case, with J± = 0.2 and Jd = 0, 
and in the frustrated case at the maximally frustrated 
point ( this point is defined below), with Ji = 0.2 and 
Jd = 0.114, for a 12 x 13 lattice. The extrapolated value 
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m/Jd 


0. 


0.114 


32 


0.43796 


0.42898 


64 


0.43928 


0.42917 


96 


0.43970 


0.42920 


oo 


0.44051 


0.42921 


QMC 


0.44075 





TABLE I: Ground-state energies as function of m and Jd for 
Jx = 0.2. 



of Eq in the unfrustrated case is in good agreement with 
the QMC. It can be seen that Eq converges faster in the 
frustrated case. The strategy in this study is to increase 
J± while the ratio Jd/J± remains close to 0.5. We do 
not expect a significant decrease in accuracy for larger 
J± because as we will see below, the interchain corre- 
lations decay exponentially in this region. Thus despite 
its perturbative nature, the TSDMRG is able to treat 
model(^l well beyond the weak coupling regime as long 
as one stays in the vicinity of the QCP. 

The high accuracy enjoyed at half-filling is somewhat 
reduced in the doped case. The reason behind this re- 
duction is simply that the size of the total Hilbert space 
increases as one goes from two degrees of freedom per 
site to three. Since we keep roughly the same number 
of states in the two cases, the accuracy will be smaller 
for doped systems. This can be seen, for instance, when 
ni = 96 states are kept in the second step. The energy 
width of the retained states drops from AE w 1.6 at 
half-filling to AE « 0.8 when the doping is x = 0.125 for 
a L = 16 system. Another difficulty which arises upon 
doping is that one needs to target, during the first step 
of the algorithm, at least three charge sectors in order 
to keep matrix elements which involve interchain hop- 
ping. A typical situation is L = 16 and x = 0.125, where 
the corresponding ground state has 14 electrons. For 
this charge sector, five spin sectors with Sz = 0,±1,±2 
were targeted. For charge sectors with 13 and 15 elec- 
trons, four spin sectors with Sz = ±0.5, ±1.5 were also 
targeted. Unlike the half-filled case, ED was not done 
for all doped systems in the first step. For instance for 
L = 16, we kept 243 states. This number of states cor- 
responds to an exact diagonalization of 12 sites which 
means three DMRG iterations were necessary to reach 
the desired size. Hence, the truncation error at the end 
of the first step increases from zero at half-filling to about 
1 x 10~ 6 in the doped case. 

Removing or adding one electron on a finite system 
can significantly affect the density if the system size is 
not large enough. For instance, during the study of a 
L = 16 lattice with the hole density x = 0.125, the nom- 
inal ground state has 14 electrons, whereas the ground 
states with 13 and 15 electrons correspond to the hole 
densitites 0.1875 and 0.0625 respectively. Thus, adding 
or removing holes on small systems induces large vari- 
ations of the hole density, resulting in large differences 
between the ground states of the different charge sectors. 





X 


0.125 


0.25 


0.375 




J 


1.4 


1.14 


0.75 


E(N) 


-E(N + 1) 


0.1012 


0.0982 


0.1355 


E(N) 


-E(N-1) 


0.1250 


0.1658 


0.1615 



TABLE II: Hole densities, corresponding exchange couplings 
and energy differences between charge sectors for L = 16. 

These variations are the largest near half-filling for small 
values of J. As a consequence, the electrons cannot jump 
between the chains and the charge degrees of freedom 
remain decoupled for small values of the interchain hop- 
ping. It is necessary to fine tune J so that these energy 
differences are smaller or of the same magnitude as t±. 
For tj_ = 0.2, which will be used for all doped systems, 
J was chosen so that the difference 6 C between ground 
state energy of charge sectors differing by one electron 
is such that S c < 0.16. Fixing S c amounts to fixing the 
chemical potential instead of the coupling. The value of 
J which satisfies this criterion depends on the hole den- 
sity x. These values are listed in Table [H] for L = 16 
systems. 



III. FINITE SIZE ANALYSIS AT THE 
CRITICAL POINT AT HALF-FILLING 

The half-filled case corresponds to weakly coupled 
Heisenberg chains that we studied in Refi^. In that study 
it was found that the model displays two phases: a Neel 
Q = (n,n) when Jd/J± ^ 0.5 and a Neel Q — (n,0) 
when Jd/J± ^ 0.5. These two phases are separated by 
a critical point where the long-distance behavior of the 
spin-spin correlations along the chains in the 2D system 
is identical to those in the disconnected chains. In this 
study, we will show that this behavior found for small 
Jj_ < 0.2 extends well beyond the weak interchain cou- 
pling regime (up to J± = 0.6). The nature of the ground 
state in this regime, which was previously called "nearly 
independent chains" , will be further clarified by the anal- 
ysis of the interchain spin-spin correlations. Before start- 
ing the finite size analysis, it is important to stress the 
fact that the open boundary conditions (OBC) which 
were applied in most of this study create some difficulty 
in the interpretation of data. While OBC have the advan- 
tage of yielding higher accuracy in DMRG simulations, 
they artificially break the lattice translational symmetry. 
The OBC thus introduce a spurious dimerization d in 
the ID chain, as we will see below, which converges very 
slowly toward zero. To our knowledge, the exact behav- 
ior of this dimerization as a function of the system size 
has not been reported in the literature. On top of this 
difficulty, one must account for logarithmic corrections 
to the finite size quantities. For instance, fitting d with 
1/L, one finds that d extrapolates to a finite value in the 
limit L — > oo. This deviation from 0, which is the true 
limit, is the signature of logarithmic corrections. Log- 



J± 


L = 8 


L = 10 


L = 12 


L = 16 


0.2 


0.124 


0.119 


0.114 


0.112 


0.4 


0.252 


0.243 


0.236 


0.226 


0.6 


0.386 


0.372 


0.361 


0.342 



TABLE III: Values of Jd corresponding to the maximum of 
the ground-state energy for different J± and L. 



arithmic corrections in an open Heisenberg chain have 
been studied in Ref^&. They are due to marginally irrel- 
evant operators and their mathematical form is hard to 
find exactly for small systems. A functional form exists 
only in the limit of large L. For instance, the finite size 
spin gap has the following form 



A, 



nv 4irg(L) 



(9) 



where v = tt/2 is the spin velocity and g(L) a function 
such that 



9(L) 



V3 



4"7T log L 



(10) 



when L is very large. For small chains, g(L) has to be 
introduced as an unknown parameter which is obtained 
by fitting the data to known limits. Since the main in- 
terest of this study lies in the relative behavior of the 2D 
system with respect to an isolated ID chain, we will not 
use g{L) during the finite size analysis. When performing 
the finite size analysis, simple linear or quadratic func- 
tions of 1/L will be used even though they do not lead to 
the exact values in the thermodynamic limit for the ID 
system because of the logarithmic corrections. 



A. Ground-state energies 

The ground-state energies shown in Fig.([TJ show sim- 
ilar behavior for all values of Jj_ when Jd is varied. 
Starting from Jd = where the ground state is a Neel 
Q = (7r, 7r) phase, Eq typically increases with Jd until 
it reaches a maximum, labeled the maximally frustrated 
point J™ aa: . Then Eq decreases when Jd > J d nax where 
the ground state is in a Neel Q = (ir, 0) phase. The ra- 
tio J d nax /J± depends on both J± and L. For small Jj_, 
J™ ax /J^ remains very close to 0.5 but deviates from this 
value as J± increases. Previous studies ( see Refi^ for 
a review) found J d nax /J± —0.6 when J± = 1. This de- 
viation is larger for OBC than for PBC. This is because 
the maximally frustrated point corresponds to the point 
where the transverse component of the Fourier transform 
of the coupling 



J(Q) = 2 cos Q x + 2cosQ y [J± + 2JdCosQ x ] (11) 
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FIG. 1: Ground-state energies for L = 16 and Jj_ = 0.2(a), 
J± = 0.4 (b) and, J± — 0.6 (c) as function of Jd- 



vanishes when Q x — n. At this point 2Jd bonds cancel 
one J± bond. But when OBC are used on a L x [L + 1) 
lattice, one has L 2 bonds but only 2L(L — 1) Jd bonds. 
Thus for small lattices, there is a deficit in Jd bonds which 
translates into the shift of the maximum towards larger 
Jd for a fixed J±. 

This shift has been argued to be indicative of the exis- 
tence of an intermediate phase, possibly dimerized, in the 
region between Jd = 0.4 (where the Neel Q = (ir, ir) order 
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FIG. 2: Maxima of the ground-state energies as function of 
the lattice size for J± = 0.2 (circles), J± — 0.4 (squares) and, 
J± — 0.6 (diamonds). 
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FIG. 3: Ground-state energies for L — 16 and J z ± = 0.2 in 
the Ising case. 



parameter was found to vanish) and Jd = 0.6 (where the 
Neel order parameter Q = (w, 0) or (0, tt) vanishes). Our 
results shown in Table (^Qj, which list the position of the 
maxima for various Jj_ and L, indicate that this devia- 
tion of J d nax / J± from 0.5 is simply a finite size effect. For 
a fixed L, J™ ax / J ± increases with J± in agreement with 
the finding in the isotropic case. But one can see that 
for a fixed J±, J™ ax / 'J j_ decreases with L. Fig.(|5J dis- 
plays the extrapolated values of J d nax /J± for J± = 0.2, 
0.4, and 0.6. They all converge to the vicinity of 0.5. 
Aside from numerical uncertainties, it is thus likely that 
J™ ax /J± = 0.5 for all J± in the thermodynamic limit. 

The curves of Eq in Fig. ( ^) show that Eq is differ- 
entiable, which suggests that the transition could be of 
second order. This is to be constrasted with the pure 
Ising equivalent of Hamiltonian ( ^) obtained by setting 
all XY terms to 0. Eq for the Ising case shown in Fig. ( 01 
is not differentiable at J d nax /J±, as expected since this 
transition is of first order. It is to be noted that in this 
case where quantum fluctuations are absent, J d 1lax /J± 
does not occur at 0.5. This is another fact which favors 
the conclusion that the deviation from 0.5 seen in previ- 
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FIG. 4: Ci for L = 16 and J ± = 0.2(a), J ± = 0.4 (b) and, 
J± = 0.6(c) as function of ,J d . 

ous studies 24 is a finite size effect rather than some subtle 
quantum effect due to the existence of an intermediate 
phase between the two Neel ordered states. 



B. Short-distance spin-spin correlations 

The first neighbor transverse correlation function 

(Fig. ( a) 
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FIG. 5: Minima of Ci as function of the lattice size for J± 
0.2 (circles), J± = 0.4 (squares) and, J± = 0.6 (diamonds). 
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TABLE IV: Values of Jd corresponding to the minimum of 
Ci for different J± and L. 



0.001 



0.0008 



0.0006 



0.0004 



0.0002 - 




(12) 



taken in the middle of the lattice (il) — (L/2,L/2 + 1) 
is linear in the vicinity of J d nax /J± and vanishes at the 
minimum point j9. Table ( ITV|) lists the position of J° 
for various J± and L. These values of J° are equal to or 
near J™ ax found for Eq for all values of J± studied. We 
believe that the small differences between the values of 
J° and J™ oa: are due to numerical errors on C\ and Eg- 
The difference between these positions are in all cases less 
than 0.1%. Since the DMRG technique tends to have a 
better accuracy on the ground-state energies than on the 
correlation functions, the position of the maximum in Eq 
is taken as the reference for the transition. The position 
of the minimum shown in Fig. ( |SJ extrapolate as for Eg 
towards Jd/J± = 0.5. 

Since C\ = at the transition point, the spins in dif- 
ferent chains remain uncorrelated as in the classical case. 
This is indeed the best way to minimize frustration. In 
order to avoid the cost in energy induced by frustrated 
bonds, the system forms the largest unfrustrated clus- 
ter allowable. In this case, the largest unfrustrated clus- 
ters are independent chains. This state can actually be 
seen as the generalization of the Majumdar-Gosh state 
known in ID systems 23 . However, this does not mean 
that the chains are completely disconnected as in the 
classical case. In the results on Eg above we were un- 
able to find a point in the parameter space where the 
ground-state energy of the coupled 2D system is exactly 
equal to the ground-state energy of disconnected chains. 




FIG. 6: C 2 for L = 16 and J ± = 0.2(a), J± = 0.4 (b) and, 
J± = 0.6 (c) as function of J d . 



The ground state energies of the 2D systems are always 
lower than that of completely disconnected chains. This 
means that despite the fact that C\ — 0, all other trans- 
verse correlations are not necessarily zero at Jd — J° as 
in the classical case. 

The second neighbor correlation 

C'2 = (S35f I+3 ) (13) 

shown in Fig.(EJ), taken at (il) = (L/2, L/2 + 2), does not 
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FIG. 7: Value of C2 at the minimum as function of the lattice 
size for J± — 0.6 (circles). 
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FIG. 8: Evolution of the dimerization as function of the lat- 
tice size for J± — 0. (circles), J± = 0.2 (squares), J± — 0.4 
(diamonds) and, J± — 0.6 (triangles). 



vanish. It instead has a minimum at J™" for all values 
of J± ■ The values of C2 corresponding to Jd = J$ shown 
in Fig. (UJ for J± — 0.6 extrapolate to a finite value in the 
thermodynamic limit. Thus, at the point where C± = 0, 
C2 has a small finite value. It is worth noting that the 
dimerization is not the only way to avoid frustration as 
it is commonly believed^. 



C. Dimerization 

The examination of the short-range correlations above 
proves that at the maximally frustrated point, there is 
a weak bond between the chains. The ground-state in 
this region of parameters has been predicted via field 
theoretic approachesi&ii to be a columnar dimer state, 
with the dimerization occuring along the chains. In order 
to see if dimerization occurs along the chains, we compute 
the quantity 



(^2il^2i+ll) 



2i+2U 



(14) 



One should note that since we use the OBC, we are in 
the most favorable situation for the dimerization. i.e., the 
OBC breaks the translational symmetry of the chains and 
even isolated finite chains are dimerized. This spurious 
dimerization decays very slowly as L — > 00. It is thus 
expected that if indeed the system dimerizes at the QCP, 
the finite size behavior of du of the coupled system will 
be different from that of isolated chains. In Fig.([8J) d is 
shown for J± — 0, 0.2, 0.4, and 0.6 and for various sizes 
at J™ ax in the relevant cases, d converges to the same 
limit whether J± = or not. Since the single chain is not 
dimerized, this shows that the 2D system at the critical 
point is not dimerized either. 

We are aware of a different conclusion reached by a 
recent bosonization approach coupled to a renormaliza- 
tion group and mean-field analyses^. In that study for 
a two-leg ladder in the vicinity of the maximally frus- 
trated point, it was found that although the action of Jj_ 
is largely cancelled by Jd, at the second order there re- 
mains residual couplings J\ = 2J^/tt 2 and J2 = — SJ^/ir 2 
which are relevant. J\ favors Neel order and J2, the ring- 
exchange term, favors a dimerized state. These authors 
then carried out a self-consistent mean-field approxima- 
tion of the competition between J\ and J2 in the 2D sys- 
tem. They found a narrow dimerized state between the 
two ordered magnetic states. According to their conclu- 
sions, our prediction of the spin liquid state at the tran- 
sition is due to finite size effects. But our results above, 
performed for larger interchain couplings for which the 
dimerization should be more apparent and is not seen, 
rather support our conclusions in Ref— . There is addi- 
tional evidence, discussed below, obtained from the study 
of the three-leg ladder— which also agrees with the ab- 
sence a dimer phase. It is important mentioning that 
this type of RG analysis is not without a certain bias. 
The RG transformation usually leads to a proliferation 
of couplings. There is often some degree of arbitrariness 
in the choice of the relevant couplings. 



D. Long-distance spin-spin correlations 

Near the maximally frustrated point, the long distance 
behavior of the in-chain spin-spin correlation functions is 
nearly identical to that of the disconnected chains. These 
correlations shown in Fig. (01 for J± = 0.2 and L = 16 
are nearly independent of m. This is because, for a purely 
ID system of this size they are already known with very 
good accuracy. Since the chains are nearly disconnected, 
no significant increase of m is required to describe them 
with the same level of accuracy as in ID. 

This conclusion is supported by the behavior of the 
transverse correlation functions shown in Fig.( I10|l for 
the same set of parameters. C±(i) is found to be equal 
to for i > 4 for all the values of m. 

Cu(i) and C±(i) are independent of J± at the maxi- 
mally frustrated point, i.e. only the ratio Jd/ J± is impor- 
tant up to intermediate values. This is seen in Fig.( ITT|) 
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FIG. 9: Longitudinal spin-spin correlation at J d nax as func- 
tion of the distance for L — 16 and, m = 32 (circles), m = 64 
(squares) and, m — 96 (diamonds). 
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FIG. 11: Longitudinal spin-spin correlation at J™ ax as func- 
tion of the distance for J± — 0. (circles), J± = 0.2 (squares), 
J± = 0.4 (diamonds) and, Jj_ = 0.6 (triangles). 
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FIG. 10: Transverse spin-spin correlation at J™ ax as function 
of the distance for L = 16 and, m — 32 (circles), m = 64 
(squares) and, m = 96 (diamonds). 



FIG. 12: Transverse spin-spin correlation at J™ ax as function 
of the distance for J± = 0.2 (circles), J± — 0.4 (squares) and, 
J± = 0.6 (diamonds). 



and Fig.( I12fl where these quantities are compared for 
J± = 0.2, 0.4 and 0.6. Cii (i) is nearly identical to that of 
independent chains while C±(i) remains very small and 
decays exponentially with i. This behavior is valid from 
weak to intermediate coupling regimes. It is to be noted 
that since the maximally frustrated point is chosen at the 
maximum J™ ax of Eg, C\ is not equal to because the 
two quantities differ slightly as shown in Tabic ( ITTT1 II Vj 
because of numerical uncertainties. 

This behavior of C±(i) at Jd ~ J±/2 is to be con- 
trasted to that of the magnetic regime. In Fig.f I13|) we 
see that for instance, Cj_(4) of the magnetic case is al- 
ready four orders of magnitude larger than in the dis- 
ordered case. The exponential decay of C±(i) seen in 
Fig.f I13|) and the decay of Cn(i) which is nearly that of 
decoupled chains suggest that for Jd ps Ji/2, the sys- 
tem displays a spin analog of the sliding Luttinger liquid 
(SLL) discussed in fermionic models2iiS*ii. In the re- 
gion Jd ~ J_l/2 the two competing magnetic fluctuations 
(•7T, 7t) and (-7T, 0) cancel each other leading to irrelevant 



interchain couplings. These irrelevant couplings renor- 
malize the energy but do not modify the behavior of the 
correlation functions. Our findings show unambiguously 
how the SLL concept is linked to quantum criticality. It 
also provides evidence of fractionalization at the critical 
point as first suggested in Re£i£. 

E. Finite size spin gap 

Logarithmic corrections and the spurious dimerization 
introduced by the OBC is accompanied by a spurious gap 
in the pure ID system. In Fig. ( I14|) we see that at the 
maximally frustrated point, the gap for J± = 0.2 con- 
verges to the same value as that of Jj_ = 0. It is thus 
reasonable, knowing that the single chain is gapless, to 
conclude that the 2D system is also gapless in the thermo- 
dynamic limit. We note that the logarithmic corrections 
are not taken into account in the gaps predicted by exact 
diagonalization studies^. 
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FIG. 13: Transverse spin-spin correlation as lunction of the 
distance for J± = 0.4, J d = (circles), J± = 0.4, J d = 0.219 
(squares). 
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FIG. 14: Finite size gap with OBC as function of the lattice 
size for J± — (circles) and, Jj_ = 0.2 (squares). 



F. Implications for the isotropic Ji J2 model 

For the isotropic case (Jj_ = J) it is usually said that 
there is a large amount of evidence supporting the fact 
that the magnetic order vanishes for 0.4 < Jd < 0.6, i.e., 
a spin gapped state exists in this region^. But in our 
opinion, this conclusion is not that compelling. The re- 
sults displayed above indicate that the limit Jd/J± — * 0.5 
is in fact analogous to the problem of coupled unfrus- 
trated chains (Jd — in Hamiltonian( ^) in the limit 
J± — > 0. For this latter problem there was a prediction 
of a spin liquid state for J± < 0.1 from ED, field theory 
mapping and spin wave approaches^. However, a QMC 
studji2& later showed that the Neel order extends down to 
J± = 0.01. It is now believed that long-range order sets 
in as soon as Ji 7^ 0. This suggests that the disordered 
phase often claimed to exists for the isotropic J\ — J2 
model, which was obtained from ED of only a maximum 
of 6 x 6 systems or series expansions such as dimers ex- 
pansions which may be biased toward dimerization, is 
not a settled issue. 

In fact there is a strong evidence, unfortunately ne- 



glected, which points to the contrary. A DMRG study 27 
performed on a three-leg ladder found that the system 
is gapless for all values of J± and Jd- The lenght of the 
chains were up to L = 360 so that finite size effects were 
negligible. The computed phase diagram has two states 
a symmetric doublet phase and the quartet phase. The 
low energy properties of the symmetric doublet phase are 
similar to those of a single spin-1/2 chain. The scaling 
behavior of the quartet phase is identical to that of a 
spin-3/2. These two states will naturally evolve towards 
the Neel Q — (n,ir) and Q = (tt, 0), respectively, when 
the number of chains goes to infinity. This is fully con- 
sistent with our conclusions. 

The perturbative method applied in this study can- 
not be used directly to study the isotropic case, since 
starting with the decoupled chains is too biased. But if 
the principles uncovered in the anisotropic case are to be 
applied, it appears that even in the isotropic situation, 
the nearly decoupled chains ground state is still the best 
way to minimize frustration. From this consideration, we 
applied the method to the isotropic case and arrived at 
the same conclusions as above. Namely that the ground 
state at the maximally frustrated point is made of nearly 
disconnected chains and is thus a SLL. 

Another indication for this is the finding by a recent 
ED studjiSi that even at J± = 0.8 J the behavior observed 
at lower coupling still persists. However, this study also 
predicts a spin gap. We believe this apparent contradic- 
tion is due to the fact that the ED study was done with 
an even number of chains and the logarithmic corrections 
were not considered. 



IV. DOPED SYSTEMS 

Having identified the ground state at half-filling as a 
spin version of a SLL, we now turn to the doped case. 
We will mostly restrict ourselves to the vicinity of the 
maximally frustrated point Jd ~ J±/2. This is primarily 
because the results at half-filling indicate that far from 
this point no significantly new physics is likely to oc- 
cur. One would expect that the two Neel ground states 
found at half-filling away from the the QCP will natu- 
rally evolve into spin density wave (SDW) ground states 
upon doping. The really physically interesting question 
is: What is the effect of the doping upon the spin SLL? 

The main finding in the doped case is that the mag- 
netic properties of the system do not significantly change 
as long as one remains close to half-filling. But if the 
hole doping becomes too large, x w 0.3, the model dis- 
plays a qualitative change. The frustration ceases, and 
Jd cooperates instead of competing with J±. 

A. Ground-state energies 

Eq at x = 0.167 for L = 12 as function of Jd is shown 
in Fig. ( I15fl . It displays a maximum as for the half-filled 
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FIG. 15: Ground state energy for x = 0.167, L = 12, and FIG. 17: Ci for x = 0.167, L = 12, and J± = 0.4 as function 
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FIG. 16: Ground-state energy for x = 0.4, L — 10, and J± 
0.2 as function of Jd 



FIG. 18: d for x = 0.4, L = 10, and J ± = 0.2 as function of 
Jd 



case. One can note that this maximum is now shifted 
towards higher values of Jd- This maximum goes from 
jrnax = q.236 for x = to J^™ = 0.30 for x = 0.167. By 
inserting the holes, the system becomes less frustrated. 
It is therefore necessary to increase Jd in order to balance 
the action of J±. We find that if the hole density is in- 
creased above x ~ 0.3 the maximum in Eq is suppressed. 
For x > 0.3, Eg always decreases with Jd- This is seen 
for instance in Fig.f HB)l where Eq for x = 0.4 is shown. 
This shows that for large dopings, contrary to the near 
half-filling situation, Jd increases instead of reducing the 
stability of the system. 

B. Short-distance spin-spin correlations 

This difference between the near half-filling and larger 
hole densities is also displayed in C\ . In Fig. ( IT7|) corre- 
sponding to x = 0.167, C\ goes to at J" = 0.3. One 
may note that there seems to be a jump at the tran- 
sition point. This jump was not observed at half-filling. 
This could indicate a possible first order transition. How- 



ever this is not supported by the behavior of Eq shown 
above. We notice during the simulations that, if one is 
too close to the transition point, for the same set of pa- 
rameters, one can go from one side of the transition when 
the number of chains is small to the other side when the 
number of chains gets larger. We thus believe that the 
jump observed in C\ is somehow related to the reduced 
accuracy in the doped case induced by the mixture of the 
two competing ground states near the transition point. 
The regularity of the curve far from the transition point 
justify this assumption. In agreement with the behavior 
of Eg, the transition is suppressed when the hole dop- 
ing is too large. C\ shown in Fig.f ITB|) always remain 
antiferromagnetic and increases with Jd- 

C. Long-distance correlations 

C±(i) shown in Fig. f 119ft for x — 0.125 retains essen- 
tially the same behavior as in the half-filled case. It de- 
cays exponentially at the maximally frustrated point for 
J± = 0.2, 0.4, and 0.6 while the decay is significantly 
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FIG. 19: Transverse spin-spin correlation as lunction of the 
distance for J± — 0.2, Jd — (circles), and at the maximally 
frustrated point at J± = 0.2 (squares), J± = 0.4 (triangles), 
J± = 0.6 (diamonds). 
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FIG. 20: Transverse spin-spin correlation as function of the 
distance for J± — 0.2,Jd = 0.11 for x = 0.125 (circles), x = 
0.25 (squares) and x = 0.375 for L = 16. 



slower in the unfrustrated case. In the absence of frus- 
tration, Cj_(i) has an antiferromagnetic signature while 
Cii(i) is SDW like with kp = ~(1 — x), where x is the 
density of holes. Though we have not made a finite size 
analysis, it is reasonable to believe that the system is in 
an SDW ground state in the thermodynamic limit. At 
the maximally frustrated point, the magnetic order pa- 
rameter vanishes. 

As x is further increased to 0.167, 0.25, and 0.375, we 
find that the role of frustration is reduced. Ultimately, 
the system evolves towards a conventional SDW ground 
state. This is seen in Fig.( I20[l where the exponential 
decay of C±(i) is significantly slowed at large hole dop- 
ings. There is already four orders of magnitudes between 
Cx(4) for x — 0.125 and x — 0.375. It appears that the 
system tends to choose hole configurations that minimize 
frustration in order to avoid the cost in energy induced 
by frustration. It is in the vicinity of x = 0.125 that frus- 
tration is most effective, i.e., the magnetic order which 
exists in the non-frustrated case can easily be cancelled 
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FIG. 21: Critical values of Jd as function of doping. 

by frustration. Beyond this point, frustration becomes 
less and less effective. At quarter filling the system can 
choose a hole configuration so that frustration becomes 
ineffective. Such possible configurations are shown in 
Fig. (|24Jl . but we have not made a systematic analysis 
of them. Since the introduction of holes usually tends 
to reduce the order parameter, the reemergence of mag- 
netism at the maximally frustrated point by introduction 
of holes can be viewed as the manifestation of Villain's 
order from disorder effect. 

These arguments are illustrated in Fig.( l21f> where the 
value of J" is shown as function of x. J®(x) first slowly 
decreases from half-filling until about x = 0.125 and then 
sharply rises. Thus there is a point, for x between 0.25 
and 0.5, where the critical point is suppressed. 

D. Equal-time Green's function 

The transition point is also signaled in the equal-time 
transverse Green's function 

G±(i) = {c jla c] l+ia ). (15) 

Gi = G±(l) is shown in Fig.( I22|) decays from Jd = 
to Jd — J^ lax . This decay of G\ is consistent with the 
tendency to localize the particles within the chains in- 
duced by Jd- However, in the vicinity of J™ aa: , the sit- 
uation becomes less clear. G\ strongly oscillates. These 
oscillations cease when Jd is beyond the critical region 
where G\ decays again. It appears that G\ behaves dif- 
ferently in the two ordered states; it increases from the 
transition point in the first and decreases from the transi- 
tion point in the second state. The curious behavior seen 
in the transition region could well be the manifestation 
of the mixture of the two ground states seen in Eq and in 
C\ . This mixture affects the value of fermion correlation 
functions more severely than their spin counterpart. For 
this reason, we have not calculated charge density or su- 
perconducting correlations which involve four fermions. 
There is however no reason to suspect a charge density or 
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FIG. 22: Transverse equal-time first neighbors Green's func- 
tion at x — 0.167 as function of Jd for L — 12. 



O.Olr 




FIG. 23: Transverse equal-time Green's function at x = 0.125 
for Jd — 0.2, Jd = for spin up (circles), spins down (squares); 
J± — 0.2, Jd = 0.11 for spins up (triangles), spins up (dia- 
monds) 



superconducting ordering in the vicinity of J 7 d nax . These 
transition would be of first order and would thus be ap- 
parent in the ground state energy. 

Despite the uncertainties caused by the mixture of the 
two ground states on the two sides of the transition, we 
nevertheless show in Fig. ( I23fl the equal-time transverse 
Green's function in the SDW case and in the magneti- 
cally disordered case. The decay in the disordered state 
is slower than in the magnetic case. This raises the pos- 




FIG. 24: Holes (light color) configurations that minimize frus- 
tration for x = 0.25. The left illustration represents the half- 
filled case where the system is fully frustrated. 



sibility of the relevance of the single particle hopping. 
If true, this would mean that the spin and the charge 
behave differently in the transverse direction at the tran- 
sition point. More work needs to be done in order to 
clarify this point. 

V. CONCLUSIONS 

In this study, we have presented numerical evidence for 
the connection between quantum phase transition and 
Luttinger liquid physics from the study of coupled t — J 
chains. The ordered Q = ((1 — x)ir, ir) Neel state that 
exists at half-filling and near half-filling for J± ^ and 
Jd = is destroyed at J d — J™ a2 \ At this point, the in- 
chain spin-spin correlation of the 2D system decay nearly 
like those of independent chains. A careful examination 
of the transverse spin-spin correlations reveals that the 
chains are loosely bound and that these correlations de- 
cay exponentially. We thus identify the state of the sys- 
tem at this maximally frustrated point as a sliding Lut- 
tinger liquid. These properties are observed even when 
the interchain coupling are well beyond the perturbative 
regime 

The key mechanism behind this connection is mag- 
netic frustration. The leading effect which brings a com- 
promise between the frustrated bonds is the severing of 
these bonds to form the largest unfrustrated clusters with 
the lowest energy. These unfrustrated units are then cou- 
pled by some residual interactions generated by quantum 
fluctuations. These residual interactions may or may not 
drive the system away from the physics displayed by the 
largest unfrustrated clusters. This picture is in constrast 
to the widely accepted hypothesis which suggests that the 
leading mechanism to avoid frustration is short or long- 
range dimerization which may (spin-Peierls) or may not 
(Resonating Valence Bond) lead to the breaking of the 
lattice translational symmetry. It is however to be noted 
that the two pictures do not necessarily always conflict 
with each other. In some models, the largest unfrus- 
trated clusters could just be dimers or plaquettes. Our 
point here is that this is not the case in the anisotropic 
Ji — J2 model. 

In the model studied in this paper, LL physics arises at 
the QCP because the severing of frustrated bonds leads 
to nearly independent chains which are the largest un- 
frustrated clusters. These chains are coupled, even when 
the original interchain couplings are large, only by resid- 
ual interactions which do not appear to destabilize the 
ID physics at the QCP. In a numerical simulation, it is 
however impossible to rule out completely the emergence 
of relevant interchain interactions at a very low energies. 
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